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ABSTRACT 

Using data obtained by the EUV Imaging Spectrometer (EIS) onboard Hinode, we have per- 
formed a survey of obvious and persistent (without significant damping) Doppler shift oscillations 
in the corona. We have found mainly two types of oscillations from February to April in 2007. 
One type is found at loop footpoint regions, with a dominant period around 10 minutes. They are 
characterized by coherent behavior of all line parameters (line intensity, Doppler shift, line width 
and profile asymmetry) , apparent blue shift and blueward asymmetry throughout almost the en- 
tire duration. Such oscillations are likely to be signatures of quasi-periodic upflows (small-scale 
jets, or coronal counterpart of type-II spicules), which may play an important role in the supply 
of mass and energy to the hot corona. The other type of oscillation is usually associated with the 
upper part of loops. They are most clearly seen in the Doppler shift of coronal lines with forma- 
tion temperatures between one and two million degrees. The global wavelets of these oscillations 
usually peak sharply around a period in the range of 3-6 minutes. No obvious profile asymmetry 
is found and the variation of the line width is typically very small. The intensity variation is 
often less than 2%. These oscillations are more likely to be signatures of kink/Alfven waves 
rather than flows. In a few cases there seems to be a it/ 2 phase shift between the intensity and 
Doppler shift oscillations, which may suggest the presence of slow mode standing waves according 
to wave theories. However, we demonstrate that such a phase shift could also be produced by 
loops moving into and out of a spatial pixel as a result of Alfvenic oscillations. In this scenario, 
the intensity oscillations associated with Alfvenic waves are caused by loop displacement rather 
than density change. These coronal waves may be used to investigate properties of the coronal 
plasma and magnetic field. 

Subject headings: Sun: corona — Sun: oscillations — line: profiles — waves — solar wind 
1. Introduction 
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In the past 15 years there was a rapid growth of 
literature reporting quasi-periodic oscillations in 
the upper solar atmosphere (for latest reviews see, 



e.gjBaneriee et alJ2QQ7HDe Moortel 



- , . _ Nakariakov 

2012). These oscillations are often explained as 
signatures of various mod es of MHD waves. For 



example. I Jess et aTl ( 20091 ) reported an oscillation 
associated with a large magnetic bright point and 
interpreted it as torsional Alfven wave. Intensity 
oscillations of a flare loop observed in radio wave- 
lengths were interpreted as fast sausage mode and 



1 



applied to constrain the density contrast between 



the loop and the background ([Nakariakov et al. 
2QQ3t lAschwanden et al.ll2QQ4h . Evidence of fast 



important for diag n osing t h e stellar coronae (e.g 



Mitr a-Kraev et al 
20091 ). 



2QQ5C iPandev fc S rivastava 



sausag e modes was also rec ently found in magnetic 
pores ([Morton et all 1201 if ). 

Lateral transverse oscillations are often inter- 
preted as fast-mode kink waves or Alfven waves. 
The ubiquitous presence of transverse waves has 
been observed in chromqspher i c spicules and mot- 
tles (iDe Pontieu et all 1200/1 [Singh et al.l [2007 : 



201 lh and coronal loops (e.g., iBerghmans 



Zaaarashvili et al. 1200/1 Kukhianidze et al. 2006: 

Kim et al.l2QQ8HHe et alJ2QQ9bHOkamoto fc De Pontieul _ _ _ tl 

201 lj Kuridze et al. 20 13), coronal loops ([As chwandenel^ BhOll: [Marsh et al.l I2003U2009I: Iking et al.l l2QQ3t 
1999 1 : Nakaria kov et al.ll999i jNakariakov fc Ofmanl r " - TT _._..._ TT , i t, 



20011: IVerwichte et alJ 12004 iTomczvk et al. 200/1 
Tomczvk fc Mclntoshll2009l : iMcIntosh et al.llioTll 
Singh et al.l 1201 if ), and plu me- like structures in 
different regions o f the Sun ([Mcintosh et al.ll201ll : 
Tian et al.ll2011bl ). Case studies suggest that some 
transverse waves are likely associated with im- 
pulsive phenomena such as flare activities, fila- 
ment eruptions, or c orona l mass ej ections (CME) 



(lAschwanden et"al 
2QQ2t IWang et al 



1999 



2UU2: IWang et al.l 12004 lU'Bhea et al.l laiJUYbl 
Chen et al.l 120081: lOfman fc Wand 120081: iHe et al 



2002; Schrii ver et al 



P'Shea et al I l2007bt 



2009 a: Verwichte et al. 20 10U lAschwanden fc Schriiver 
201 it IPietarila et al.l lioili iLiu et aD l201lL l2Q12fc 
Shen fc Li u 2012; White & Verwichtdl2012U White et al. 



20121 : IWang et al.l l2012bl ). Large-scale trans- 
verse waves with clear signatures of damping 
have also been reported in coronal streamers 
and found to be driven by CME impingement 
([Chen et alJlioiol l201lt iFeng et al.ll201ll ). Mod- 
eling efforts have also been made to understand 
the generation and propagation of these trans- 
verse waves (e.g., iMatsumoto fc Shibatal l2010l 



Murawski fc Musielakl l2010at ISoler et al.l l201lb 



Selwa fc O fman 2010 Sel waet al.l 2011a,b) and 



we refer to a recent review of lOfma n (2009) . Ob- 
servations of these Alfvenic waves ([Goossens et al.l 
2009), either persistent or impulsively excited, 
have been widely used for the derivation of coro- 
nal Alfven speed and magnetic field. For reviews 
of the observations of these transverse waves and 
their a pplications in coronal se is molog y, we re- 
fer to iNakariakov fc Verwichtel (120051) . lErdelvi 



and Zaaaras hvili et al 



In spite 

of the implications of such transverse waves and 
their role in the MHD seismology of the solar 
corona, the detection of such waves may also be 



Longitudinal oscillations have also been inten- 
sively studied. Using mainly broadband imag- 
ing observations, quasi-periodic upward propa- 
gating disturbances with a period of three to 
thirty minutes we re frequently obs e rved i n po- 
lar plumes (e.g., lOfman et al.l fl997i Il999l |2Qq5 
DeForest fc Gurmanl 1 19981 : Krishna Prasad et al 



Clettd 

1999: IDe Moortel et al.ll2000L 120021: iRobbrecht et al.l 



De Pontie u 2003; Ugarte-Urra et a l. 2004a. bllDe Pontieul 



2005; McEwan & De Moortel 2006; Tian & Xia 



2008 : Srivastava et al.l 20081: Srivastava fc Dwivedi 



20101: IStenborg et all l201lL lYuan fc Nakariakov 



20121 ). Propagating disturbance with a s horter pe- 



riod (5 s) has been recently reported bv lLiu etaL 



in a dark cavity region. Longitudinal dis 
turbances have also been investigate d through 
spect r oscopic observations (e.g., 



Baneriee et al 

20001: lO'Shea et al.l l2007al iBaneriee et al.l 120091: 



Gupta et al.l2QQ9U2Qlft IWang et al.ll2012aHKrishna Prasad et al 
2Q12h . which often reveal a correlation between 



changes of the line intensity and Doppler shift 



(Brvnild sen et al.lll99 9: Wikst0l et al. 2QQ ftlWang et al 
2009aUbl:lKitagawa et al.ll2010l : lMariska fc Muglach 
2010l : iGupta et al.l l2012h . These long-lasting 
propagating disturbances often have a speed of 
50-200 km s -1 and are usually interpreted as 
slow-mode magneto-acoustic waves propagating 
into the co r ona a long the magnetic field lines. 
De Pontieul ( 20041 ) proposed that some of these 



longitudinal oscillations might be caused by the 
leakage of p-modes into higher layers of the so- 
lar atmosphere along inclined magnetic field lines. 
Longitudinal oscillations with clear s ignatures of 
damp i ng hav e also b een detected ([Wang et al.l 
120021 . l2003al lbL 120051 : iMariskaet al.l 120081 ). In 
some cases there is a tt/2 phase shift between 
the oscillations of th e intensity and Dopple r shift 
([Wang et al.l l2003al lbl: iMariska et all l2008f ). sug- 
gesting the possible presence of standing slow 
waves. There are suggestions that these longi- 
tudinal oscil lations can be used to infer thermal 
conduction (Ofman & Wang 2002), adiabatic in- 
dex dVan Doorsselaere et al.l 120111 ) and magnetic 
field ( IWang et al.l 120071 ) in coronal loops. For 
reviews of the observations and theories of these 
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long-lasting longitudinal propagating disturbances 
and qui ckly damped long it udinal oscillati o ns, we 
refer tolDe Moortell (EoQ 9h iTarovan et al.l (|2QQ9h , 
iBaneriee et all <|201lh and lWand (|201lh . 

However be Pontieu fc Mcintosh] (|2010h and 



Tian et al.l (|2Qlla| ) argued that quasi-periodic sig- 



nals are not necessarily waves. They found that 
in loop footpoint regions all of the four line pa- 
rameters (line intensity, Doppler shift, line width 
and profile asymmetry) show coherent behav- 
ior, which can be easily explained by the sce- 
nario of recurring upflows (or quasi-periodically 
enhanced and wea kened upflows). In addition, 
Tian et al. (|2Qllah found a net blue shift and 



blueward asymmetry of the line profiles at the 
boundary of an active region (AR) throughout the 
entire observation duration, favoring the interpre- 
tation of a recurring high-speed (50-150 km s _1 ) 
outflow superimposed on the nearly static coro- 
nal background. The presence of these up- 
flows has been confirmed thro ugh detailed anal- 
ysis of emission line p r ofiles (|Hara et al. 20081: 

r\_ T^i j_* _j_ _ l I lorvrv/ni K J _T__ j l_ 0_ T"\ _ T- I 



The discovery of these rapid quasi-periodic up- 
flows (or jets, plasma ejections) from spectro- 
scopic observations challenges the universal wave 
interpretation for coronal oscillations, specifically 
for some longitudinal propagating disturbances. 
There is no doubt that both waves and flows ex- 
ist in the upper solar atmosphere. But it is dif- 
ficult to distinguish between them without an- 
alyzing all line parameters including intensity, 
Doppler shift, line width and profile asymmetr y 
(|De Pontieu fc McIntoshll2Qld : iTian et al.lkoilat ). 
although realistic numerical MHD mode l s ma y 
help resolve this ambiguity (| Of man et al. 2012 ). 
In this paper, we perform a survey of obvious and 
persistent (without significant damping) Doppler 
shift oscillations in the corona by using the data 
obtained by the EUV Imaging Spectrometer (EIS, 
Culhane et al. 1 l2007h onboard Hinode. Through 



detailed analysis of the temporal evolution of all 
the line parameters, we conclude that some oscil- 
lations are better explained by recurring upflows 
and others are more likely to be real waves. 



De Pontieu et al.l 120091 iMcIntosh fc De Pontieul 

2009allr]:[Peterll2010l:lBrvans et al.llioTcilMartmez-Sykora eg all Data selection and reduction 

Doha fc Zhu kov 2Qlll: IUgarte-Urra fc Warrenl 



2011 



2011 



ITian et al. 2QTT3Q, l2012t IPoschekl l2Q12h . 



Tian et aD ( 2012 ) also reported similar high- 



speed upflows in CME- induced dimming regions. 
Based on the similarities of the upward propagat- 
ing disturbance s in plu mes and AR boundaries , 
Mcintosh etail (|2010at ) and ITian et al.l (j2011bh 



suggested that the propagating disturbances in 
plumes might also be dominated by flows, al- 
though an accurate investigation of the line pro- 
file asymmetry in coronal holes is impossible with 
current instruments. If this scenario is correct, 
these high-speed episodic upflows/jets might be 
an efficient means to provide heated mass into 



De Pontieu et al. 20091: 


2009b 


; lHansteen et al.l 


2011; 


Tian et al,l l2011bh 



Mcintosh fc De Pontieu 



tions and X-Ray jets have also b een identified 



from recent imaging observations (Morto n et al. 



20121 : IZhang et al1l2Q12L ISu et al.ll2012l ). Several 



numerical studies have been performed to inves- 
tigate the formation of these plasma eject a (e.g. 



Murawski & Zaaarashvili 2010b; Murawski et al 



2011 



2011 



EIS made a large number of sit-and-stare (fixed 
slit location) observations from February to April 
in 2007. We examined the quicklook Doppler- 
grams (automatically generated level 2 data) of 
the Fe XII 195. 12 A line, which are available on the 
website of Hinode Science Data Center Europe, for 
these observations. In Table [T] we list the starting 
time of all observations in which we found obvi- 
ous (easily identified by eye) and persistent (with- 
out significant damping) Doppler shift oscillations. 
The locations of oscillations on the slits are also 
listed in Table [TJ Typical cadence of these obser- 
vations is 30 s. 

After standard correction and calibration of the 
EIS data, a running average over 3 pixels along 
the slit was applied to the spectra to improve 
the signal to noise ratio. The spatial offset in 
the Solar- Y direction between the two CCDs have 
been corrected using the standard EIS routine 
eis-ccd-offset.pro. This routine measures the off- 
set by co- aligning images from Fe Vlll 185.21 A 
and Si vn 275. 35A. We mainly selected three 
strong emission lines in each observation for our 



Srivastava fc Murawski 2011: Mar tmez-Svkora^D stud y : Fe XII 195.12A, Fe XIII 202.04A, and 



McLaughlin et al.ll201 



Fe xiv 264. 78A. The latter two lines are be- 
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lieved to be clean (jYoung et al.l 120071 ), although 
an unidentified line might be present in the red 
wing of Fe XIII 2Q2.04A ([Mcintosh et al.l l2010ri 
iTian et all l2Q12h . The Fe XII 195.12A line is 
known to be blended with the Fe XII 195. 18A 
line that i s usually a few percent o f the 195. 12 A 



line tnat is usuall y a tew percent o i tne Wo.iZJ\ 
inten sity ([Young et al.l 120071 120091 : iBrown et all 



20081 ). We used this line since it is the strongest 
EIS coronal line in non-flare conditions. The 
Fe XII 195. 12 A line was also involved in density 
diagnostics and the coalignment between TRACE 
and EIS. More lines were selected for a compre- 
hensive analysis of the 2007 March 28 observation. 

As a common practice, we applied a single 
Gaussian fit to each EIS spectrum to derive the 
line intensity, Doppler shift and line width. Since 
there are no cool chromospheric lines in the EIS 
spectral range, we simply assumed a zero average 
Doppler shift of each line in each observation. In 
some observations there was signi ficant jitter in 
the so lar-y direction. Similar to iKitagawa et al. 
(2010), we performed a cross-correlation between 
intensities of Fe XII 195. 12 A at different expo- 
sures to evaluate and remove the jitter in the 
y-direction. It turned out that the jitters were 
mostly smaller than 2" . Since the oscillations 
we found are associated with specific locations on 
the slit, we can safely conclude that they are not 
caused by possible jitters in the x-direction. We 
also calculated the average RB asymmetry in the 
velocity range of 60-120 km s _1 for each line pro- 
file. The RB asymmetry techniqu e was first intro- 
duced b y iDe Pontieu et al.l ([20091 ) and later mod- 
ified bv ITian et all (|2011ch . The RB asymmetry 
here refers to RBp, where the spectral position 
corresponding to the maximum intensity is used a s 
the line centroid (see details in lTian et al.ll2011ch . 



3. 2007 March 28 oscillations: two types 
of oscillations observed on the same slit 

Starting from 18:37 on 2007 March 28, EIS per- 
formed sit-and-stare observation in an AR contin- 
uously for more than three hours. Figure G^A) 
shows an image of this AR in TRACE 195A pass- 
band and the location of the EIS slit as determined 
by cross-correlating the EIS Fe XII 195. 12 A line 
intensity along the slit and the TRACE intensity 
at different x- locations. From the temporal evolu- 
tion of the Doppler shift of Fe xm 202. 04A shown 



in Figure DJB), we can clearly see two locations 
(indicated by the two arrows) with obvious oscil- 
lations throughout the entire duration. One oscil- 
lation (around y=162") is characterized by some 
inclined elongated features of enhanced blue shift. 
This oscillation appears to be associated with the 
footpoints (or lower parts) of AR loops, as can be 
seen from Figure QJ A). We call it oscillation 1 in 
the following discussion. The other oscillation (ar- 
round y=-95") is dominated by small red shift and 
we can clearly see the quasi-periodic appearance 
of vertical features with enhanced red shift. From 
Figure [TJ A) we can see that this oscillation (oscil- 
lation 2) appears to be associated with the tops 
(or upper parts) of a bundle of AR loops. The slit 
is found to be almost parallel to the loop plane 
and cover at least 45" of the loop top. Note that 
although the edge the bright region around y=- 
70" is caught by the slit, this part does not show 
significant oscillation and it is not analyzed in the 
following. 

We have performed plasma diagnostics for these 
two regions revealing clear Doppler shift oscilla- 
tions. As an example, Figure QJC)&(D) show the 
electron densities and differential emission mea- 
sure (DEM) curves at the positions of y=162" and 
y=-95 // . The densities were derived from the in- 
tensity ratio of the line pair Fe xn 196.64/195. 12A. 
The theoretical relationship between the line ra- 
tio and d ensity, as extracted from the CHIAN TI 
database (Dere et al. 1997; Landi et al. 2006) 



is 

shown as the black line. The electron densities at 
the two locations are log (A^/cm _3 )=9.54±0.12 
and 9.20±0.04, respectively. More lines were used 
for the derivation of DEM curves. We averaged the 
profiles of each line acquired at different exposures 
and assu med an uncertainly of 10% for the in- 
tensities flTrip athi et al. 2008). By using the rou- 
tine chiantLdem .pro (also used by lLee et al. I l201ll 



Tian et alllioiih available in SolarSoft (SSW) and 



assuming a constant pressure of 10 cm K, we 
obtained the DEM curves at the two locations. 
The peak temperatures at the two locations are 
around log (T/K)=5.85 and 6.15, respectively. 

3.1. Oscillation 1: correlated changes of all 
line parameters 

These oscillations are also found in several other 
emission lines with different formation temper- 
atures. From the left panels of Figure [2] we 
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(A) TRACE 1 95 intensity (arbitrary unit) ( B ) Doppler shift of Fe XIII 202.04 (km/s) 

: i , : : 

0.00 0.17 0.33 0.50 0.67 0.83 1.00 -30 -20 -10 10 20 30 




200 250 300 350 400 450 500 50 100 150 

X (arcsec) Minutes after 1 8:37 on 2007/03/28 




0.02 0.04 0.06 0.08 0.10 0.12 4.5 5.0 5.5 6.0 6.5 

Ratio of Fe XII 196.64 / (195.12+195.18) Log (T/K) 



Fig. 1 — (A) A TRACE 195 A image taken at 20:34 on 2007 March 28. The black vertical line indicates the 
location of the EIS slit. The arrows point to the approximate locations of two persistent oscillations. (B) 
Temporal evolution of the Doppler shift of Fe xill 202. 04 A along the slit. The short horizontal bars in (A) 
and dashed lines in (B) are used to indicate the ranges of illustration in Figs. 2-5. (C) Electron densities at 
the positions of y=162 // and y=-95 // , respectively. Error bars indicate the standard deviations of the density 
time series. (D) DEM curves at the positions of y=162 // and y=-95 // , respectively. 
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Fig. 2. — Left: Temporal evolution of the Doppler shift of several emission lines in the range of y=130 // ~210" 
in Figure QjB). Right: The black, red, green, and violet curves represent respectively the detrended line 
intensity, Doppler shift (inverted), line width, and RB asymmetry (inverted) averaged over the region between 
the two dashed lines shown in the left panels. The intensity has been normalized to the local background 
and is shown as the percentage divided by 2. The RB asymmetry has been normalized to the peak intensity 
of the line profile and is shown in the unit of percentage. The units of Doppler shift and line width are 
km s _1 . For the purpose of illustration, the intensity and RB asymmetry are offset by 3 and -4 respectively 
on the y-axis. 
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can see that oscillation 1 is clearly present in 
all lines with a formation temperature of 0.6-2.0 
MK. The inclined elongated features of enhanced 
blue shift are most clearly seen in the strongest 
line Fe XII 195. 12 A and can also be identified in 
other weaker lines. For all coronal lines, we see 
blue shift throughout almost the entire observa- 
tion duration. The transition region (TR) line 
Si VII 275. 35 A reveals a mixture of both blue 
shift and red shift, which might be related to 
the enh ancement of downflows at lower temper - 
atures (|Tian et al.l l2Qlld : iMdntosh etHI l2Q12h . 
We averaged all the line parameters (intensity, 
Doppler shift, line width and RB asymmetry) 
over the 11 pixels between the two dashed lines 
and produced a time series of each line param- 
eter. By subtracting a 10-mi nute running aver - 
age from each time series (e.g. JWang et al.ll2009af ) 
we obtained a time series of each detrended line 
parameter, which is shown in the right panels 
of Figure [2j We found that the intensity vari- 
ation generally follows the variation of the line 
width. It is also clear that the Doppler shift 
and RB asymmetry show coherent behaviors. By 
taking the inverted values of Doppler shift and 
RB asymmetry (now positive values mean blue 
shifts or blueward asymmetries), we found that 
all of the four parameters show coherent varia- 
tions. Note that the less-than-ideal correlations 
at some instances are probably partly caused by 
the poor spectral resolut ion and photon noise 
of the EIS instrument ( De Pontieu fc Mcintosh! 

20ld iTian et al.l 1201 lah . As demonstrated by 

^ 1 1 noise 



De Pontieu fc Mcintosh! (|2010h . the photon 



has progressively worse effects as one goes to 
higher moments, i.e., the asymmetry suffers much 
more from noise than the intensity. The fore- 
ground and background emission in the line of 
sight (LOS) direction might be another reason for 
the reduced correlations. 



can clearly see obvious blueward asymmetry in the 
loop footpoint region. A comparison between pan- 
els (D) and (F) suggests that the Fe xm 202.04A 
line is possibly affected by an unidentified blend 



r (De Pontieu fc M 


clntoshl 


2010ri ITian et al. 


2012). 



This effect tends to shift the RB asymmetry to 
the redward and thus reduce the magnitude of 
blueward asymmetry in the loop footpoint re- 
gion. However, we still used the Fe xm 202. 04 A 
line for our detailed analysis since it is much 
stronger than Fe xiv 264. 78A. The similar be- 
havior of different coronal lines as shown in Fig- 
ure [2] suggests that the effect of the possible 
blend or background issue is not significant in 
our study. Panel (L) presents a typical line pro- 
file of Fe xiv 264. 78A (nine profiles averaged) and 
we can clearly see an enhancement of the blue 
wing. Such an asymmetric line profile suggests at 
least two emission components filar a et al. 2008; 



20ld:lMcIntosh_fc De Pont ieu 2009a.b; Pet erl2Ql 



De Ponti eu et al.l 200 



np onents (Mara et al. zUUoc 
it iDe Pontieu fc Mclntoshl 
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mm 



Brvans et al.l 120101: i Martme z-Svkora et al.l 120111 
Doha fc Zhii^2QllUUgarte-Urra fc WarrenllioTll 
Tian et al.l l2011al ld: lDoschekl l2Q12) and a study of 
center to limb variation clearly suggests that the 
enhancement in the blue wing i s not caused by 
blends or noise ( Tian et al 1 l2012h . The correlated 
changes of different line parameters can be seen 
from the detrended time series in both panels (E) 
and (G). Panel (J) and (K) show the cross correla- 
tions between the Doppler shift and intensity/line 
width/RB asymmetry as a function of time lag. 
The nearly zero time lag for each of the parame- 
ter pair confirms the coh erent behavi o rs of a ll line 
parameters. Similar to I Wang et al. ( 2009a ). we 
measure the oscillation period as the value cor- 
responding to the peak global wavelet power and 
the uncertainty as the half width at half max- 
imum (HWHM). The oscillation periods of the 
intensity and Doppler shift are 9.01±2.64 min and 
10.7±2.70 min respectively. The Doppler shift 
period of 10.7±2.70 min is consistent with the 
typical time gap between two inclined stripes in 



The detailed temporal evolutions of all of the 
four line parameters of the strong Fe xm 202. 04 A 
line are presented in Figure (3^A)-(D). The loop 
footpoint region is clearly characterized by quasi- 
periodic enhancement of the line intensity, blue 

shift a nd lin e width, sim ilar to what be Pontieu fc Mclntostfe left panels of Figured The large uncertainties 
(|2010h and ITian et al.l ([2011ah found. We can of periods suggest that the time gap between two 

also see weak signatures of blueward asymmetry events (i.e., upflows, see below) varies significantly 

in the region between the two dashed lines. The with time. The oscillation amplitude, defined by 

Fe xiv 264.78A line is clean and from panel (F) we the square root of the peak global wavelet power 
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Fig. 3. — (A)-(D) Temporal evolution of the line parameters (intensity: detrended, others: original) of 
Fe xm 202.04A in the range of y=130 // ~210 // in Figure [T^B). (E) Temporal evolution of the detrended line 
intensity (i), Doppler shift (v), line width (w), and RB asymmetry (a) averaged over the region between the 
two dashed lines shown in panels (A)-(D). The line styles, units, and illustration methods are the same as any 
right panel in Figure O The root-mean-square values of these parameters are also shown to the right. (F) 
Similar to (D) but for Fe xiv 264.78A. (G) Similar to (E) but for Fe xiv 264.78A. (H)-(I) Wavelet spectra for 
the detrended intensity and Doppler shift of Fe xm 202. 04 A. The periods and amplitudes derived from the 
global wavelets are indicated as Pi/Pv and Ai/Av, respectively. The dashed lines indicate a significance level 
of 99%. (J)-(K) Cross correlations between the Doppler shift and intensity(black)/line width(green)/RB 
asymmetry (violet) for Fe xm 202. 04 A and Fe xiv 264. 78 A. The maximum correlation coefficients (CC) 
and the corresponding time lags (shift) are also marked. (L) An observed Fe xiv 264. 78 A line profile and 
measurement errors are shown as the diamonds and error bars, respectively. The green line is the single 
Gaussian fit. The two dashed red lines represent the two Gaussian components and the solid red line is 
the sum of the two. The velocity (v) and exponential width (w) derived from the single (SGF) and double 
(lst/2nd for the two components) Gaussian fits are shown in the panel. Also shown is the intensity ratio of 
the two components (a). 



(|Wang et al.ll2009al) , is 6.25% for the intensity and 
1.67 km s _1 for the Doppler shift. Note that the 
data gaps separate the entire time series into three 
segments and that we selected the longest segment 
for wavelet analysis. We have also performed the 
wavelet analysis for the other two shorter segments 
and the obtained periods are close to that of the 
central segment. 

3.2. Oscillation 2: oscillation most clearly 
seen in Doppler shift 

Oscillation 2 can be clearly identified in all lines 
with a formation temperature ranging from 1.3 
to 2.0 MK. From the left panels of Figure |4] we 
can see that oscillation 2 is characterized by al- 
most vertical stripes of enhanced red shift. To 
the best of our knowledge, such obvious long- 
lasting oscillations in warm (1.3-2.0 MK) coronal 
lines have never been found in observations before. 
Lower-temperature lines such as Si VII 275. 35 A 
and Fe x 184. 54 A are weak in this loop top re- 
gion so that we could not study the detailed tem- 
poral evolution of their Doppler shift. Similar to 
Figure [2j we obtained time series of the detrended 
line parameters for each line in the region between 
the two dashed lines and present them in the right 
panels of Figure HI It is clear that the oscillation is 
most prominent in Doppler shift. The variations 
of the intensity, line width and RB asymmetry are 
all typically smaller than those of oscillation 1. 

The temporal evolution of all of the four line 
parameters of Fe xill 202. 04 A is detailed in Fig- 
ure GJA)- (D). Clearly, oscillation 2 is mainly 
present in the Doppler shift. No clear discernible 
variations are found in the line width and RB 
asymmetry. The Fe xiv 264. 78 A line shows be- 
haviors similar to Fe xill 202. 04 A, as can be seen 
from panels (F) and (G). The line profiles in this 
loop top region are typically symmetric and an ex- 
ample is presented in panel (L). The variation of 
the detrended intensity is smaller than that of os- 
cillation 1. But from panel (G) we see a seemingly 
correspondence between the intensity and Doppler 
shift variations. The correspondence can also be 
seen from Figure |4] and is more clear in hotter 
lines such as Fe xiv 274. 20 A and Fe XV 284. 16 A. 
From the analysis of cross correlation, we find a 
maximum correlation coefficient at the time lag 
of about 0.9 min between the time series of in- 
tensity and Doppler shift for the two lines shown 



in Figure [U Such a time lag corresponds to a 
phase shift of about tt/2 or 1/4 period, given 
a Doppler shift oscillation period of 4.91±0.67 
min. However, since the correlation coefficient is 
smaller than 0.3 (which can easily be caused by 
random, chance correlations) and the dominant 
period (12.7 mini 1.01) of the intensity is largely 
different from that of the Doppler shift, we can 
not make any solid conclusion about it. The os- 
cillation amplitude is 1.78% for the intensity and 
1.77 km s _1 for the Doppler shift. Clearly, the 
intensity fluctuation is much smaller than that of 
oscillation 1. Whilst the amplitude of the Doppler 
shift oscillation is very similar to that of oscillation 
1. 

4. Statistical results and discussion 

Table [T] lists the observational information and 
analysis results, including observation time, lo- 
cations on the slit, root-mean-square values of 
all line parameters, oscillation amplitudes, peri- 
ods and uncertainties of the intensity and Doppler 
shift, electron densities, and locations, for all os- 
cillation events we identified. These parameters 
were derived from an averaged time series over 
11 pixels along the slit in each event. We care- 
fully selected these 11 Y-pixels in each event and 
made sure that 1) the temporal behavior at these 
11 pixels is almost the same, 2) this 11-pixel re- 
gion covers the major oscillating part. Similarly, 
we detrended each time series by subtracting a 
10-minute running average. We mainly used the 
strong Fe xill 202. 04 A line which is available in 
all observations for our statistical studies. The 
weaker but clean Fe xiv 264. 78 A line was also 
analyzed for a confirmation of the characteris- 
tics revealed by Fe xill 202. 04A. The line pair of 
Fe xii 196.64/195. 12A or Fe xn 186.88/195. 12A 
was used to diagnose the electron densities. Im- 
ages obtained by the Tra nsition Region and Coro- 
nal Explorer (TRACE, lHandv et all Il999h , the 
X-Ray Telescope (XRT. iGolub et al.l 120071 ) and 
the Extreme Ultraviolet Im aging Telescope (EIT, 
Delaboudiniere et al.| [l995) were used for identify- 
ing the coronal structures associating with these 
oscillations. Based on these statistics we can sub- 
divide the oscillations we found here into two 
types. One type (type I) are mainly found in the 
lower part of AR loops (Sect 4.1) and they are 
probably related to recurring upflows. The other 
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Table 1: Persistent Doppler shift oscillation events observed by EIS from Feb to April in 2007. The following 
information is listed for each event: observation date (yyyymmdd), starting time (hhmmss), lowest y-pixel 
of the selected 11-pixel region on the slit (Y), root-mean-square values of the detrended intensity (Ri, %), 
Doppler shift (Rv, km s _1 ), line width (Rw, km s _1 ) and RB asymmetry (Ra, %), intensity oscillation 
period and uncertainty (Pi & Piun, min), Doppler shift oscillation period and uncertainty (Pv & Pvun, 
min), amplitudes of intensity (Ai, %) and Doppler shift (Av, km s _1 ) oscillations, electron density and 
uncertainty (N & Nerr, log cm -3 ), oscillation type and location. 



id date 



time Y Ri Rv Rw Ra Pi Piun Pv Pvun Ai Av N Nerr Type, location 



20070201 
20070202 
20070203 
20070220 
20070221 
20070326 
20070327 

8 20070328 

9 20070328 

10 20070216 

11 20070221 

12 20070223 

13 20070223 

14 20070224 

15 20070224 

16 20070225 

17 20070326 

18 20070327 

19 20070328 

20 20070328 

21 20070419 



013212 351 
004912 245 
005642 193 
175013 296 
021812 310 
150012 168 
143412 220 
031531 264 
183726 303 
163920 345 
194813 394 
125943 384 
185342 380 
005742 381 
130643 370 
112113 185 
150012 29 
160926 247 
031531 12 
183726 40 
191102 189 



1.07 0.44 
0.89 0.49 
1.14 0.59 
4.57 2.79 
3.45 1.62 
0.87 0.48 
1.24 0.56 
1.84 0.89 

2.08 0.56 
1.93 0.70 
1.22 0.54 

1.09 0.46 
1.06 0.56 
1.27 0.54 
0.88 0.63 
0.91 0.66 
0.79 0.68 
0.80 0.74 
0.84 0.72 
0.64 0.67 
3.17 0.97 



0.30 0.46 
0.51 0.65 
0.56 0.63 
2.33 2.60 
1.33 3.12 
0.48 0.51 
0.41 0.49 
0.77 0.81 
0.77 0.76 
0.46 0.71 
0.34 0.56 
0.22 0.45 
0.35 0.55 
0.34 0.57 
0.25 0.31 
0.29 0.46 
0.27 0.38 
0.28 0.38 
0.27 0.46 
0.26 0.33 
0.44 0.65 



10.7 2.44 9.82 
7.57 3.19 10.7 
8.77 3.00 7.57 

11.6 2.04 12.7 

10.7 2.90 9.82 
8.26 1.88 9.82 

13.8 2.42 9.82 
13.8 3.17 11.7 
9.01 2.64 10.7 
5.35 1.29 4.13 
10.7 1.34 3.47 
6.37 1.46 4.91 
6.14 0.87 3.47 
5.84 1.68 3.18 
4.50 1.03 3.47 
7.57 3.06 4.91 
8.26 2.92 5.35 
6.94 2.05 4.91 
8.26 3.44 7.57 
12.7 1.01 4.91 
6.94 1.02 5.35 



2.05 
4.79 
2.44 
2.45 
2.05 
2.05 
4.89 
2.99 
2.70 
1.46 
0.43 
2.19 
0.61 
0.67 
0.30 
1.29 
0.93 
1.12 
2.59 
0.67 
1.29 



3.89 
2.36 
2.71 
14.0 
11.8 
2.53 
4.11 
5.50 
6.25 
4.15 
4.73 
2.09 
2.38 
2.48 
1.67 
1.95 
2.21 
2.04 
2.03 
1.78 
8.39 



1.65 9.23 
1.45 9.25 
1.60 9.22 
8.47 9.65 
5.19 9.32 
1.45 9.21 
1.40 9.65 
2.63 9.03 
1.67 9.52 
1.43 9.67 

1.57 8.83 
0.95 9.27 
1.15 9.25 
1.08 9.31 
1.76 9.32 

1.58 9.10 
2.13 9.00 
1.70 9.16 

1.76 8.94 

1.77 9.14 
2.73 9.08 



0.04 
0.04 
0.04 
0.13 
0.17 
0.05 
0.10 
0.05 
0.07 
0.07 
0.05 
0.07 
0.08 
0.09 
0.09 
0.05 
0.04 
0.03 
0.04 
0.03 
0.03 



, loop leg 
, loop leg 
, loop leg 
, loop leg 
, loop leg 
, loop leg 
, loop leg 
, loop leg 
, loop leg 
I, QS BP 
I, loop top 
I, loop top 
I, loop top 
I, loop top 
I, loop top 
I, unclear 
I, loop top 
I, loop top 
I, loop top 
I, loop top 
I, QS BP 



10 



Doppler shift (km/s) 



-60 
-80 
100 
120 
-60 
-80 
100 
120 
-60 
-80 
100 
120 
-60 
-80 
100 
120 



20 -10 
^Fe XII 195.12 



10 



20 



^Fe XIII 202.04 



=Fe XIV 274.20 



^Fe XV 284.16 



50 100 150 

Minutes after 18:37 on 2007/03/28 




Intensity (% X 0.5) 

Doppler shift (km/s) 

Line width (km/s) 

RB asymmetry [60-120 km/s] (%) 



J I I I I I I I I I I I I I 



_ ;W ■ ' 

-6 | | | | | | | | | L 

6 
4 
2 





I I I - 

I I L 



I I L 







50 100 150 

Minutes after 18:37 on 2007/03/28 



Fig. 4.— Similar to Figure[2]but for the range of y=-133"~- 53" in Figure[U;B). The Doppler shift and RB 
asymmetry values are not inverted here, which is also different from Figure [2] 



type of oscillations (type II) are found mostly near 
the top of loops, and the observational facts seem 
to suggest an interpretation of transverse wave. 
We have to mention that event 6 is very close to a 
localized bright region (TRACE 195A passband) 
in the core of an AR. But from the XRT image 
it is clear that the oscillating part (the localized 
bright region in TRACE 195A passband) is the 
leg (or lower part) of a hot loop system. We could 
not find any context images for event 16 so that 
its location is unclear. 



4.1. Type I oscillations: 
flows? 



recurring up- 



We found nine oscillation events which reveal 
properties similar to oscillation 1 discussed in 
the previous section. One example is shown in 
Figure [6f A). These oscillations are all found at 
the foot point regions (lower parts) of AR loops. 
They are characterized by coherent behaviors of 
all line parameters (line intensity, Doppler shift, 
line width and profile asymmetry), obvious blue 
shift and blueward asymmetry throughout almost 
the entire observational duration. The time lag 
between the Doppler shift and other parameters is 
very close to zero and their correlation coefficients 
are typically larger than 0.3 for all of these events. 



The line widths are obviously enhanced at the os- 
cillation locations. These oscillations are named 
type I oscillations in the following discussions. 

For these type I oscillations, the global wavelet 
power is often distributed over a wide range of 
periods. The dominant periods (period corre- 
sponding to the peak of a global wavelet) range 
from 7.57-13.8 min for the intensity and 7.57- 
12.7 min for the Doppler shift, with an average 
of 10.5±2.3 min for the intensity and 10.3=bl.4 
min for the Doppler shift. The oscillation am- 
plitudes are in the range 2.36-14.0% for the in- 
tensity and 1.40-8.47 km s _1 for the Doppler 
shift. The average amplitudes of the intensity and 
Doppler shift are 5.9% and 2.8 km s _1 , respec- 
tively. In two events star ting; at 17:50:13 on Feb 
20 (previously ana l yzed b y lTian et all ( 2011a ) and 
Nishizuka fc Haral feullh ) and 02:18:12 on Feb 21 



we see an intensity amplitude of about 13% and 
Doppler shift amplitude of about 7 km s _1 . The 
root mean square values of the line width and 
RB asymmetry for these two events are about 2 
km s _1 and 3%, respectively. These values are 
much larger than those of other events. The 
magnitudes of fluctuations are likely influenced 
by effects such as LOS projection and the fore- 
ground/background emission. 



11 



> -120 ^ 



>■ 



-80 
100 
120 

o " 60 

§ " 80 
j6 -100 

>- -120 

o " 60 

§ " 80 
js -100 

>- -120 





(C) Fe XIII 202.04 
Width (km/s) 



=(D) Fe XIII 202.04 

L RB Asymmetry [60-120 km/s] (% 



CD CD 

=1 

cdE 

£: CO 

cB 



CD 



CD 



|E) FIB XIII 202.04 



z(F) Fe XIV 264.78 

D RB Asymmetry [60-1 20 km/s] (%] 



6 
4 
2 

-2 
-4 
-6 
-60 

-80 

jg -100 

>- -120 



§|G) PeXIVl264.78 

CDCD 4 s - 7 

?o 2 

i§ .8 



Ri-0.64 
Rv=0.67 

lw=0.26 
Ra=0.33 





110 


1-10 



Ri=1 .32 
Rv=0.81 

Rw=0.37 
Ra=0.41 



50 



100 

Minutes after 18:37 on 2007/03/28 



150 



(H) Fe XIII 202.04: Wavelet of Intensity 




(I) Fe XIII 202.04: Wavelet of Doppler shift 



100 
Time (min) 

(J) Fe XIII 202.04: Cross Corr. 



14QD.0 1.6 3.2 
Power 




100 
Time (min) 



140.0 1.6 3.2 
Power 



g °- 5 

1 0.0 

CD 

S "0.5 



(K) Fe XIV 264.78: Cross Corr. 



-4 



■i shift=0.88CC=0.14 

-2 2 
Time lag [Min] 



o 



0.5 
0.0 
-0.5 




(L) Fe XIV 264.78: Line Profile 



— v-ishift=0.93CC=0.15 

-4-2 2 
Time lag [Min] 



SGF J\ 




v:1.79 / 




w:47.5 f 









264.6 264.7 264.8 264.9 265.0 
Wavelength (angstrom) 



Fig. 5. — Similar to Figure [3] but for the range of y=-133 // ~-53 // in Figure DJB). In panels (J)-(K) only the 
results of cross correlation between the Doppler shifts and intensities are shown. In panel (L) only the single 
Gaussian fit is applied to the observed line profile. 
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Fig. 6. — (A) The Doppler shift and time series of all line parameters of Fe xiv 264. 78A for a type I 
oscillation event. The line styles, units, and illustration methods are the same as in Figure 03(G). (B)-(D): 
The Doppler shift and time series of all line parameters of Fe xill 202. 04 A for three type II oscillation events. 
The line styles, units, and illustration methods are the same as in Figure E](E). (E)-(F) Wavelet (of Doppler 
shift, similar to Figure Efl)) and cross correlations analysis (similar to Figure [5fJ)) for the event shown in 
panel (D). 
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The line profiles are usually enhanced in the 
blue wings, suggesting the existence of at least 
two emission components. Compared to the sin- 
gle Ga ussian fit, the RB-guided double Gaus- 
sian fit foe Pontieu fc McIntoshll2Qld: iTian et all 



l2011al R l2012h usually does a much better job 
From Figure |3]^L) one can also conclude that 
the enhanced line width and significant blue 
shift derived from the single Gaussian fit are 
actually caused at least partly by the super- 
position of a weak high-speed upflow compo- 
nent on a strong b ackground emission component 



(|Tian et al.ll2011ch . 



These observational facts point to the expla- 
nation of recurring high-speed upflows (or quasi- 
periodically enhanced and weakened upflows) for 
type I oscillations, as p reviou sly su ggested by 



type l osculations, as p revio usly suggest ed by 
iDe Pontieu fc Mclntoshl (|2Q10h and ITian et al. 



(|2011ah . Figure [7] presents synthetic line pro- 



files and time series of the line parameters in the 
scenarios of recurring upflows and periodically en- 
hanced and weakened upflows. The background 
coronal emission component and high-speed up- 
flow component are represented by the two red 
dashed lines. Clearly, continuous high-speed up- 
flows with periodic enhancement of the flow inten- 
sity can easily explain the quasi-periodic enhance- 
ment of the line intensity, line width, blue shift 
and blue ward asymmetry. In real observations, 
the flow intensity might be enhanced differently 



in dif ferent periods (clearly seen in ITian et al 



which can sometimes cause a deviation of 
the dominant period derived from wavelet analysis 
from the real period of recurrence. 

Correlated changes between intensity and 
Doppler shift are usually interpreted as propa- 
gating slow mode ma gnetoacoustic waves (e.g ., 



Wang et al.l l2009al lbl: iNishizuka fc Haral I2Q11D . 



However, here we can see that such a correla- 
tion also exists in the case of recurring upflow. 
It seems that the line width and asymmetry are 
the ke y to distinguish b e tween slow waves and up- 
flows. Verwi chte et al. I (|2010ah argued that slow 
waves can also cause line asymmetries when the 
emission line is averaged over an oscillation pe- 
riod or when a quasi-static plasma component in 
the line of sight is included. However, the fre- 
quency doubling of the line width, as shown in 
their Figure 3, is not seen in our type I oscilla- 
tions. In addition, the irregular quasi-periodic 



changes (change differently in different periods) 
of the intensity and Doppler shift are also ex- 
pected in the scenario of recurring upflows since 
each repetitive upflow is likely to be independent 
of each other. So at present we conclude that 
the type I oscillations we identified here are more 
likely to be dominated by quasi- periodically en- 
hance d and weakened upflows . iHeggland et al. 



(|2QQ9h and iMcLaughlin et al.l (|2Q12h have 



re- 



cently demonstrated that such quasi-periodic out- 
flows can be generated by os cillatory recon n ection 



(McLaughlin et al 2009^ iMorton et al (2012) 



also proposed that the recurrent plasma ejec- 
tions they observed are reconnection-driven and 
similar to type-II spicules. Bas ed on an ob- 
servation of a supersonic blob, ISrivastava et al 



(|2Q12I ) suggested that multi-temperature plasma 
blobs could be generated by a recurrent three- 
dimensional (3-D) reconnection process via the 
separator dome below the magnetic null point, 
between the emerging flux and pre-existing field 
lines in the lower solar atmosphere. Indeed, quasi- 
periodic flow pulses associated with a null point 
have recently be en reported (jSu et al. 12012 ). On 



the other hand, iTripathi et al. ( 2012 ) explained 
the hot plasma upflows observed in the warm 
loops as evidence of chromospheric evaporation 
in quasi- static coronal loops. Whilst in th e 3-D 
simulation of Martmez-Sykora et al. ( 2011 ). the 
chromospheric plasma is ejected as result of being 
squeezed by the magnetic tension. 

We would like to point out that slow waves may 
also exist. The disappearance of enhanced line 
width and blueward asymmetry at higher parts of 
loops might be caused by LOS effect or low sig- 
nal to noise ratio of EIS . But it might also be ex - 
plained by sl o w wa ves ([Nishizuka fc Haral 120 lib . 
lOfman et al.l (|2Q12h recently found that period- 
ically driven upflows can excite undamped slow 
waves in loops. They showed that the properties 
of quasi-periodic flows are manly related to driv- 
ing source, while the properties of waves are deter- 
mined by the structures of loops. It is likely that 
high-speed upflows dominate in footpoints (lower 
corona) of some loops and slow waves dominate 
higher up in the loops. More recently, through a 
detailed analysis of the temperature dependence 
of the speeds of t he propagating disturbances, 
Kiddie et al.l (|2Q12h found signatures of both up- 
flows and slow waves. On the other hand, slow 
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Fig. 7. — Synthetic line profiles (diamonds) and time series of the line parameters in the scenarios of 
recurring upflows (A) and periodically enhanced and weakened upflows (B). In the panels of line profiles, 
the green/red solid line is the single/double Gaussian fit. The two dashed red lines represent the primary 
(background) and secondary (upflow) components, respectively. The red and green vertical lines indicate 
the line center of the primary component and single Gaussian fit, respectively. The single Gaussian fit 
parameters (i-line intensity, v-Doppler shift, w-line width) and RB asymmetry (a) are shown in each panel. 
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waves associated with the leakage o f p-modes may 



well e xist in the chromosphere (e.g. Jde Wijn et al. 



2009) and some of them might propagate further 
up into the corona and mix with the upflows. 

These high-speed upflows in lower parts of 
loops may play an important role in the supply 
of mass and energy to th e hot corona and prob- 
ably also the solar win d (|De Pontieu et al . 2009: 



De Pontieu et al.l 120111 ). There are suggestions 
that the blue shifts of t he order of 30 km s -1 at 
boundaries of some ARs (iMarsch et al]l2004ll2 008: 
Del Zann j2008UTripathi et al.ll2009UMurrav et al.l 
20101 : IWarren et all l201ll lYoung et al.l l2012h . as 
derived from a single Gaussian fit, are signatures 
of the nascent slow solar wind dSakao et al.l 120071: 
Harra et al.l 120081: iDoschek et alTl2008l : iHe et al. 



2010UBrooks fe WarrerJl201ll : ISlemzin et al.ll2Q12t 
Zangrilli fe Polettol I20I2I ). This might be true 
for the upflows along open field lines originating 
from the AR boundaries. But we have to point 
out that the speeds of these outflows are around 
100 km s _1 instead of 30 km s _1 , as revealed 
by the R B asymmetry analysis and double Gaus- 
sian fit (lHara et al.l 120081: be Pontieu et al.l 120091: 



De Pontieu fc Mcinto sh 2010: lMcIntosh fc De Pontieu 
2009a , b: Peter 2010; B rvans et alE oiO: Martmez^Sykora 
201l[|Dolla fc Zhukovll201ll:IUgarte-Urra fc Warre: 
201ll : lTian et al.ll2011aMl2012l ). In addition, some 
of these outflows are clearly associated with legs of 



actually found in quiet- Sun bright points, which 
are also believed to consist of a group of minia- 
ture loops. The oscillation is most prominent in 
Doppler shift and sometimes the cross correlation 
analysis reveals a tt/2 phase shift between the in- 
tensity and Doppler shift oscillations. The varia- 
tions of the line intensity, line width and RB asym- 
metry are all typically smaller than those of type 
I oscillations. The average root-mean-square val- 
ues of the detrended line intensity /line width/RB 
asymmetry are 1.9%/0.83 km s _1 /l.ll% for type 
I and 1.2%/0.31 km S -70.48% for type II oscilla- 
tions. 

One type II oscillation event (event 17 in Ta- 
ble [1]) presented in Figure |6fB) is characterized by 
almost vertical stripes of enhanced red shift. In 
addition, the slit is found to be almost parallel to 
the loop plane and covered at least 50" (coherence 
scale) of the upper segment of the loop. Such char- 
acteristics are very similar to the events 20 (see 
Figure [I] & Figured]) and 19 (not shown here) and 
the coherence scales are around 50". There are sig- 
natures of correspondence between the intensity 
and Doppler shift variations. The time lag be- 
tween the two suggests a possible tt/2 phase shift, 



closed field lines (man et al.l l2011d : iBoutrv et al. 



20121 ) and they are followed by cooling down- 
flows ( lUgarte-Urra & Warren] 2011 : Kamio et al 



20111 : Mcintosh et al.l 120121 ) instead of propagat- 
ing into the interplanetary space. These out- 
flows may be the coronal counterpart of type-II 
spicules which could supply he ated mass to the 
corona (jDe Pontieu et al.l 120111 ) . Taking an av- 
erage speed of 100 km s _1 and assuming a den- 
sity of log (A r e /cm _3 )=8, the mass flux density of 
these outflows is estimated to be about 3.3x10 
g cm" 
0.2. 



-10 



s 1 if we use a temporal filling factor of 



4.2. Type II oscillations: waves? 

Twelve oscillation events are found to reveal 
characteristics different from those of type I os- 
cillations. These oscillations are named type II 
oscillations in the following discussions. These os- 
cillations are found to be associated with the up- 
per parts or tops of AR loops. Two of them are 



efekhough sometimes the correlation coefficient is 
smaller than 0.3 and the dominant period of in- 
tensity is larger than that of the Doppler shift. In 
all of these events, we do not see clear signatures 
of propagating features (inclined stripes) in the 
time-slit diagrams of intensity and Doppler shift. 
Such an observational fact suggests that the phase 
speed is at least of the order of 1200 km s _1 (50" 
divided by the exposure time 30 s) if these oscilla- 
tions are signatures of propagating features. Such 
a high speed is much larger than the sound speed 
and comparable to the Alfven speed in the corona. 

The event 11 (presented in Figure EJC)) only 
exhibits oscillation in a narrow region (~10") on 
the slit, which might suggest that the loop width 
is about 10" since the slit was found to make a 
large angle with respect to the loop top (see Fig- 
ure [8j A)). No clear correlation was found for the 
intensity and Doppler shift variations in this event. 
We have found that the top parts of loops associ- 
ated with oscillating events 12, 13, 14, 15 and 18 
all made a relatively large angle with respect to 
the EIS slit. The coherence scales of these oscil- 
lations are typically 10"-30". These type II oscil- 
lations are usually associated with quasi-periodic 
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Fig. 8. — Coronal structures associated with two type-II oscillation events (events 11 & 21). (A) An XRT 
image taken at 19:49 on 2007 Feb 21. The black vertical line indicates the location of the EIS slit. The 
region between the two horizontal bars corresponds to the y-range of the Dopplergram in Figure [6fC). The 
arrow points the approximate location of oscillation event 11. (B) An XRT image taken at 17:53 on 2007 
April 19. The region between the two horizontal bars corresponds to the y-range of the Dopplergram in 
Figure [6jD). The arrow points the approximate location of oscillation event 21. 
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appearance of (enhanced) red shift and no signif- 
icant blue shift has been found, which is a puzzle 
but might be partly related to the assumption of 
zero Doppler shift of the average line profile. 

Most of these type II oscillations are found in 
ARs. However, two events (events 10 & 21) are 
found to be associated with coronal bright points 
(BPs) . In these two events we see a relatively large 
intensity fluctuation and a clear ir/2 phase shift 
between the intensity and Doppler shift oscilla- 
tions. Analysis results of event 21 are presented 
in Figure El^D) -(F) and the oscillating BP can be 
found in Figure E^B). The blue shift between y=- 
20"~5" and red shift between y=-45 // — 20" might 
be related to our assumption of zero Doppler shift 
of the average line profile. It could also be caused 
by a siphon flow along the small BP loop system. 

The periods of the type II Doppler shift oscilla- 
tions are in the range of 3-6 min, with a mean 
of 4.6 min. The global wavelets usually peak 
sharply and thus the uncertainties of the peri- 
ods are smaller compared to type I oscillations. 
The average oscillation amplitude, as computed 
from the global wavelet, is 1.6 km s _1 . Such a 
small oscillation amplitude is extremely difficult to 
be identified from on-disk imaging observations if 
they correspond to transverse displacement oscil- 
lations. With spectroscopic observations, we can 
clearly see them. 

Unlike type I oscillations, it is difficult to find 
any significant profile asymmetry at the oscilla- 
tion locations. The almost symmetric line pro- 
files seem to exclude the presence of high-speed 
outflows which are believed to be responsible for 
the type I oscillations. Perhaps a possible mech- 
anism for these type II oscilla tions is the inter- 
mittent obscuration by fibrils (|De Pontieul l2005t 
de Wijn et al.l 120071 ). or the coronal response of 
fibril motion. Imagining cool chromospheric fibrils 
going up and down quasi-periodically, the overly- 
ing hot coronal plasma can be forced upward and 
released downward with the same period. Such 
a scenario seems to be similar to vertical piston 
motion. In this case the entire coronal plasma is 
moving at a small speed, which may explain the 
absence of profile asymmetry and line broadening. 
However, it is unknown whether the fibril motion 
can affect the tops of these coronal loops which 
are often high (e.g., ~50" for the type I oscillation 
event shown in Figure [T]) in the corona. 



4-2.1. Mode identification and coronal seismol- 
ogy 

More likely, the type II oscillations are signa- 
tures of MHD waves. In the following we discuss 
possible wave modes which might explain our ob- 
servations and their roles in coronal seismology 
when applicable. 

Propagating slow wave 

Slow mode magneto-acoustic waves are compress- 
ible modes and a correlation (or anti-correlation) 
between intensity and Doppler shift is expected in 
the propagating slow wave. Our observations re- 
veal very small intensity fluctuations of all type II 
oscillation events except events 10 and 21. More- 
over, we do not see any clear signature of corre- 
lation (or anti-correlation) between intensity and 
Doppler shift oscillations in any of our type II os- 
cillation events. Thus, we can rule out the possi- 
bility of propagating slow waves. 

The top parts of oscillating loops are more or 
less parallel to the EIS slit in events 17, 19 and 
20 (see the example in Figured]). Given a typi- 
cal coronal sound speed of 200 km s _1 , the wave 
would need 180 s (much larger than the cadence 
of 30 s) to propagate through a distance of 50" 
so that we should be able to see propagating fea- 
tures (inclined stripes) in the time-slit diagrams of 
intensity and Doppler shift. The absence of such 
propagating features again rules out the possibil- 
ity of propagating slow waves. 

Standing slow wave 

The seemingly tt/2 phase shift between inten- 
sity and Doppler shift in some oscillation events 
suggests the possib il ity of standing slow waves 
([Sakurai et al.l l2QQ2t lltitagawa et all l2010h . For 
the fundamental mode the wave period is deter- 
mined by twice the loop length divided by sound 
speed ([Wang et al.ll2002[ ). Due to the strong fore- 
ground emission and the mixture with emission 
from surrounding structures, it is often difficult 
to accurately locate the legs of the oscillating 
loops (so loop reconstruct i on tec hni ques are highly 
desire d, e.g., iFeng et all ( 2007 ) & Svntelis et al. 
([2012[ )). Supposing that the oscillating loop legs 
are located at [300",-130"] and [320 // ,-30 // ] in Fig- 
ure Q] and assuming a semi-circular loop, we obtain 



18 



a loop length of 160". Using a sound speed of 200 
km s _1 , the period was estimated to be 19.22 min, 
which is much larger than the observed period 
(4.91 min) of event 20. Therefore, fundamental 
mode of standing slow waves can not explain type 
II oscillations in such large loops. A possibility 
that can be considered is a higher harmonic. But 
compared to fundamental mode, harmonic modes 
are difficult to be excited in coronal loops. In addi- 
tion, for most of the type-II oscillations the corre- 
lation between intensity and Doppler shift is poor 
(coefficient smaller than 0.3), which again does not 
favor a slow wave interpretation. 

However, standing slow waves may not be ex- 
cluded for oscillations found in smaller loops such 
as the events of Feb 16 and April 19 (events 10 
& 21 in Tabled! Figure Ep)-(F)). In these two 
events we also see relatively larger fluctuation of 
intensity, similar periods of intensity and Doppler 
shift oscillations, and a tt/2 phase shift between 
intensity and Doppler shift, which favors the in- 
terpretation of standing slow waves if the intensity 
fluctuations are signatures of density ch ange (e.g., 
ISakurai et al.ll2002l iKitagawa et al.lEoiol ) . In such 
small BP loop systems, undamped propagating 
waves from one loop leg might reflect when they 
reach the other loop leg, thu s producing standing, 
waves. The first oscillation in lErdelvi et al.l ( 20081 ) 
might be of similar type. However, as we will ex- 
plain in the following, the intensity fluctuations 
might be related to the fact that different parts of 
a loop are sampled at different times. In this case, 
the intensity fluctuations are not related to density 
change but may be caused by loop displacement 
during transverse oscillations. 

Standing kink wave 

For the fundamental mode of a standing kink wave 
the period is determined by twice t he loop length 



V Ai 



(2) 



divid e d by the kink speed (e.g., Robert s et al 
1984 : IWang et all l2QQ2t lAschwanden fc Schriiver 
201 if ) defined as: 



Ck 



■Vm 



(1) 



where pi and p e are internal and external densi- 
ties, respectively. Va% is the internal Alfven speed 
which is related to the internal magnetic field {Bj) 
and density in the form of: 



In order to explain the period (4.91 min) of oscil- 
lation event 20, the kink speed needs to be about 
775 km s -1 . This value is in the range (543-2322 
km s _1 ) of those measured from TRACE trans - 
verse oscillations ( Ofman fc Aschwandenl [2002). 
However, we have to bear in mind that the value of 
Ck depends on the magnetic field strength and den- 
sity contrast between the loop and background, 
both of which can vary a lot in different cases. 
In addition, the oscillation could be modulated by 
some continuous driver like p- modes, as suggested 
by the fact that the oscillations are persistent and 
undamped. Thus, the period may differ from the 
value predicted by twice the loop length divided 
by the kink speed. Also, the above phase speed 
is based on cylindrical or slab geometry approxi- 
mation of real coronal loops. Finally, kink modes 
can also exhibit some intensity disturbances if the 
angle between the LOS and t he oscillating struc- 
ture is not exactly 90 degree ([Cooper et al.l 12003 : 



Wang et al.l l2012bh . So the standing kink mode 



can not be ruled out. The weakly damped Doppler 
shift oscillati on event (2nd event ) interpreted as 
kink mode bv lErdelvi et al. ( 2008f ) may be of simi- 
lar type. They did not find any flare or prominence 
eruption associated with the oscillation, which is 
also similar to our type II oscillations. 

If these waves are standing kink waves, we 
can in principle obtain the internal Alfven speed 
and magnetic field strength according to Equa- 
tions [TB421 Using the line ratio method, we have 
calculated the internal and external densities from 
observations. However, the density ratio p e /pi we 
obtained is usually in the range of 0.2-0.7, which 
is much larger than the density ratio normally 
assumed (~0.08). Such results can be explained 
by the fact that LOS summ ation substantially re- 
duces the density contrast (|De Moortel fc Pascod 
20121 ). Fortunately, the density enters through 
square root and the dependence of the phase 
speed on the density and the contrast is weak. 
In fact there is only about 25% difference in 
Ck between the density ratios of 0.08 and 0.7. 
Here we simply took the typi cal value of 0.08 
(Aschwanden fc Schriiverl l201lh and obtained a 



value of 570 km s 1 for the internal Alfven speed 
in oscillation event 20. Taking the calculated den- 
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sity of log (A/ i /cm _3 )=9.14, the magnetic field 
strength of the oscillating loop is estimated to 
be 11.6 G. Here the internal mass density pi is 
calculated as 1.2m p Ni, where m p is the proton 
mass and the constant of 1.2 is due to the con- 
sideration of Helliu m abundance (e.g. JPriestll 19821 : 
Wang et al 1 l2009al ). Using a density uncertainty of 



loosely the magnetic oscillations we observed with 



log (TV^ /cm -3 ) =0.03, a period uncertainty of 0.67 
min, and assuming a 20% uncertainty of the loop 
length , the uncertainty of the est imated magnetic 
field (|Nakariakov fe OfmanllioOlh would be -25% 
(2.9 G). Note that Equation [T] is based on the 
simple configuration of a thin and straight cylin- 
drical tub e embedded in a uniform magnetic at - 
mosphere (Ed win fc Robertsl 

Eas3 [Roberts 2008), 
which might be too idealized for our coronal loop 
oscillat ions. In addition, through nu merical simu- 
lations [De^MooiteL^£ascoi ( 20091 ) have demon- 
strated that the magnetic field derived from coro- 
nal seismology might differ from the real magnetic 
field by —50% in some cases. 

Propagating Alfvenic wave 

Ubiquitous coronal Alfven waves have been found 
by observations of bot h the Coronal Mult i- channel 

m 200l 



Polarimeter (CoMP, 



Tomczyk_elaL| [2007 



iTomczvk fc McIntoshl[200i)_ and the Atmosph 



eric 



Imaging Assembly (AIA, iLemen et al.l l2Q12h on 
board the Solar Dyn ami cs Observatory (SDO) 
(Mcintosh et al. | l201lh , IVan" iDoorsselaere et al. 



20081 ) proposed that these waves are better called 
fast mode kink waves rather than Alfven waves. 
Howev er, another theoretical paper [Goossens et al. 
(|2012l ) presents a different view and clearly states 
that these waves can be considered as surface 
Alfven waves or Alfvenic waves in the nomencla- 



ture of iGoos sens et ajj (|2 



Mcintosh et al.l 



(|201ll ) also pointed out in the online supple- 
mentary information that: 1) in the highly dy- 
namic atmosphere there are no stable cylindri- 
cal w aveguides with straight magne tic field lines 

for 



that 
kink 



Van Doorsselaere et al.l ([2008) require 

Alfven wave is corn- 



waves; 2) the term 
monly used by the communities of fusion plasma 
physics and space physics for a largely incompress- 
ible transverse wave for which the major restoring 
force is the magneti c ten sion. Here we follow 
Goossens etafl (|2009h and lMcIntosh et al.l (|201lh 



EIS. Here we follow [ Goossens et all ([20091 ) and 



iMdntosh et al.l ([201ll ) and use the term Alfvenic 
waves to describe more loosely the magnetic oscil- 
lations we observed with EIS. 



.Tomczvk et all (|2007f ) and lTomczvk fc Mcintosh 
(2009) bund that these transverse waves are most 



clearly present in Doppler shift and that their os- 
cillation period peaks around 5 min, which are 
very similar to most of the type II oscillations we 
found here. Moreover, both these Alfvenic waves 
and our type II oscillations are persistent and no 
significant damping is observed. In addition, the 
presence of the almost vertical strips of Doppler 
shift features in both Figure [5] and Figure |6fB) 
indicates that the oscillations should have a very 
large propagating speed. Given a phase speed of 
the order of 1200 km s _1 , the wave only needs 30 
s to propagate through a distance of 50". This 
duration is comparable to the exposure time and 
observational cadence, which can easily explain 
the vertical strips. All of these observational facts 
suggest that some, if not all, of the type II Doppler 
shift oscillations might be spectroscopic signatures 
of the propagating Alfvenic waves. The periods 
of these oscillations, 3-6 min, suggest that the p- 
modes may greatly modul ate the generation o r 
propagation of these waves (jTomczvk et al ] l2007l ). 
The amplitudes of these type-II Doppler shift os- 
cillations are usually in the range of 1-2 km s _1 , 
which is not that far from the velocity amplitude 
found bv lMdntosh et al.l (|201lh in AR using AIA 
data (5±5 km s _1 ). Also we have to bear in mind 
that the Doppler shift oscillation amplitudes ob- 
served by EIS are likely to be reduced due to 
the LOS integration and the coarser spatial res- 
olution compared to AIA. For a quantification of 
the effect of spatial resolution on the measured 
Doppler shifts from Alfvenic waves, we refer to 



Mcintosh fc De Pontieul (|2012[ ). 



If these transverse oscillations are propagating 
Alfvenic waves, we can still calculate the inter- 
nal magnetic field if the Alfven speed is known. 
Unfortunately, the wave only needs about 30 s to 
propagate through a distance of 50" as we men- 
tioned above. Such a short time prevents us from 
deriving the phase speed from the time-distance 
diagram (jTomczvk et al.ll2008[ ). 



and use the term Alfvenic waves to describe more 
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As we discussed above, various observational facts 
seem to favor the interpretation of Alfvenic waves 
(kink or Alfven waves) for our type II oscillations, 
although slow-mode standing waves could not be 
excluded for events 10 & 21. 

The WKB energy flux of these t ransverse waves 
can be estimated as following (e.g. . iTomczvk et al 



2007 . 120081 : lOfman fc Wand 120081 : Mcintosh et al 
20111 ): 



pi(Sv) 2 V pha 



(3) 



Taking the observed value of the velocity ampli- 
tude 5v=1.77 km s _1 , the calculated number den- 
sity of log (7V^/cm _3 )=9.14, and the derived phase 
speed of 775 km s _1 , we obtain an energy flux of 
6.7 x 10 3 erg cm -2 s _1 for event 20, which is two to 
three orders of magnitude lower than the energy 
flux required to balance the radiative and conduc- 
tive losses of the active corona (2xl0 6 erg cm -2 
s _1 ). But as we mentioned above, the Doppler 
shift oscillation amplitudes are likely to be greatly 
reduced due to the coarse spatial resolution and 
LOS integration. This is supported by recent 
observation showing close st rands in a loop sys - 
tem oscillating out of phase (|Wang et al.ll2012bh . 
Recent high-resolution observations of SDO/AIA 
clearly reveal an o scillation amplitude o f the or- 



der of 20 km s _1 (Mc intosh et al.l 1201 if ). Such a 



large amplitude leads to an energy flux compara- 
ble to that required for heating the quiet corona 
and solar wind. 

4.2.2. Intensity oscillations associated with Alfvenic 
waves 

Weak intensity oscillations could also occur in 
the case of Alfvenic waves, since periodic loop dis- 
placement could lead to the scenario that differ- 
ent parts (with different intensity) of a loop are 
sampled periodically. If the polarization of the 
swaying motion is such that both a LOS compo- 
nent (visible in Doppler shifts) and a component 
in the plane of the sky occur, then loops moving 
into and out of a spatial pixel could lead to slight 
intensity fluctuations that are not a signature of 
density changes. 

Such a scenario can be easily modeled with a 
toy model. Figure [9] shows how the intensity fluc- 



tuation, as well as a tt/2 phase shift between in- 
tensity and Doppler shift can be produced. We as- 
sume a loop carrying a simple harmonic Alfvenic 
wave with an amplitude of 3 km s _1 and a pe- 
riod of 5 min (panel C). Such a wave will lead to a 
swaying motion which causes a periodic transverse 
displacement of the loop by up to ~287 km (panel 
B). As a result of this displacement, different parts 
of the loop with different intensities are sampled 
by the instrument slit at different time (panel A). 
Periodic displacement leads to periodic sampling 
of the same part of the loop, thus producing pe- 
riodic variation of the observed intensity. We de- 
rived intensity profiles across different loops from 
several raster scans of EIS (observed on 2007 Feb 
21, 2007 May 19 and 2007 Aug 23) and present in 
panel (D) a typical normalized intensity profile of 
a loop observed by EIS. Using this intensity pro- 
file and assuming an angle of 60° between the po- 
larization direction (X in panel A) and LOS, we 
obtain an intensity fluctuation of ~1.5 % (panel 
E) if the targeted pixel is ~2" away from the loop 
center. The Doppler shift oscillation has an ampli- 
tude of ~ 1.5 km s _1 and reveals a clear tt/2 phase 
shift with respect to the intensity oscillation. 

Intensity oscillations are usually believed to be 
signatures of compressible waves based on the 
assumption that the intensity c hange is a re- 
flecti on of density change (e.g., iKitagawa et al.l 
l201Qh . A 7r/2 phase shift between intensity and 
Doppler shift is often believe d to be an indicator 
of standing slow wav es (e.g., ISakurai et al. 20021 : 
IKitagawa et al.ll2010f ). However, from our Figure [9] 
it is clear that intensity fluctuations might also be 
caused by sampling changing portions of a loop, 
with different densities, due to the periodic trans- 
verse motion. In this case, a tt/2 phase shift be- 
tween intensity and Doppler shift is expected if 
the parts of loop sampled by the instrument slit 
have a monotonically decreasing or increasing in- 
tensity profile. Note that the tt/2 phase shift here 
is caused by the fact that different parts of the 
loop are sampled at different time, and is certainly 
not an intrinsic property of the transverse wave. 
We emphasize that for mode identification when 
using ElS-like spectroscopic observations the tt/2 
phase shift alone is not sufficient and one needs 
to know information such as the associated loop 
structure and phase speed (see section 4.2.1) for a 
more definite identification. 
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Fig. 9. — Intensity oscillations associated with Alfvenic waves. (A) Transverse displacement of a loop. 
X and Y are directions across and along the loop, respectively. The region between the two dashed lines 
(-0.5"~0.5") indicate the slit location (or targeted pixel, 1" wide). The slit is assumed to be parallel to 
the loop. The dotted line marks the location of the loop center. The color coding is shown in panel (D). 
(B) Loop displacement relative to X=2". The period is 5 min. (C) Loop velocity. (D) Intensity profile 
across the loop. The part of the loop between the two dashed lines indicates the slit location at time=1.25 
min. (E) Time series of the intensity fluctuation and Doppler shift when the angle between the polarization 
direction and LOS is 60°. The time lag between intensity and Doppler shift oscillations and their maximum 
correlation coefficient are also marked in the panel. An online movie (ml.mpg) is associated with this figure. 
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The amplitude of intensity oscillation depends 
on the slope of the intensity profile in the parts 
of the loop sampled by the instrument slit and 
the angle between the polarization direction and 
LOS. The online movie (ml.mpg) shows how these 
two factors impact the amplitudes of intensity 
and velocity oscillations. Using a fixed angle be- 
tween polarization and LOS (60°), the first part 
of the movie shows that the intensity fluctuation 
becomes larger when the sampled part of the loop 
has a steeper intensity gradient. Note that the 
red line marks the central location of the sampled 
parts of the loop. If this sampled location is the 
loop center, the intensity fluctuation is small and 
the period of the intensity oscillation (if signifi- 
cant) is half of the velocity period, since the in- 
tensity slightly decreases on both sides of the loop 
center. The second part of the movie clearly shows 
that the intensity amplitude increases as the an- 
gle between polarization and LOS increases from 
0° to 90°. The velocity amplitude changes in the 
opposite sense. Except for the cases of 0° and 90°, 
the tt/2 phase shift is always there. 

As we mentioned above, the intensity oscilla- 
tion and the tt/2 phase shift between intensity 
and Doppler shift variation are much more promi- 
nent in hotter lines such as Fe xiv 2 74. 20 A and 
Fe XV 284. 16 A (e.g., Figure [p. This might be 
explained by considering the possible different in- 
tensity profiles across the loop for different lines. If 
the average temperature of the surrounding coro- 
nal plasma is comparable to the formation tem- 
perature of Fe XII 195. 12A or Fe xm 202.04A, the 
intensity profile should be flatter for these cooler 
lines and thus the possible intensity fluctuation 
associated with the transverse oscillations will be 
smaller. For hotter lines, there will be a larger con- 
trast (steeper intensity profile) between the loop 
and the surrounding (cooler) corona and hence the 
intensity fluctuation will be larger. 

The maximum correlation coefficients are al- 
ways 0.98 in our idealized toy model. In real ob- 
servations, the sampled parts of the loop may not 
always have a smooth monotonically decreasing or 
increasing intensity profile. This is likely to be the 
case as a loop is believed to consist of several sub- 
resolution strands. This effect will de-smooth the 
intensity time series and thus reduce the correla- 
tion between the intensity and Doppler shift. In 
addition, we assume that the oscillating part of 



the loop is parallel to the instrument slit, which 
is certainly not always the case in observations. 
If the loop makes a large angle with respect to 
the slit and there is no change of density along 
the loop, the observed intensity change will be re- 
duced and the tt/2 phase shift might be less clear. 
In fact we see a more obvious it/ 2 phase shift in 
events 10, 17, 19, 20 and 21. In all of these events 
the top part of the oscillating loop is more or less 
parallel to the EIS slit. In other type-II oscilla- 
tion events the oscillating loop makes a large an- 
gle with respect to the slit and we generally do not 
see a clear tt/2 phase shift. Another major factor 
that would cause a reduced correlation, a variable 
phase shift and intensity oscillation amplitude is 
the possible change of the sampled loop location 
throughout the timeseries. This could happen due 
to either natural evolution of the loop or slight 
wobble/jitter in the EIS pointing. Finally, the 
instrument noise will certainly reduce the corre- 
lation coefficient and cause a variable phase shift 
with time. 

5. Conclusion 

We have performed a statistical study of 
Doppler shift oscillations using data taken by HIN- 
ODE/EIS. We have found mainly two types of os- 
cillations: one type (type I) is mainly found at loop 
footpoint regions and the other (type II) is typi- 
cally associated with the upper part of loops. The 
type I oscillations generally show coherent behav- 
ior of all line parameters (line intensity, Doppler 
shift, line width and profile asymmetry), apparent 
blue shift and blueward asymmetry throughout 
almost the entire observational duration. These 
observational facts seem to be consistent with 
the scenario of episodic high-speed upflows (fine- 
scale recurrent jets). The type II oscillations, 
with a period of 3-6 minutes, are most clearly 
seen in the Doppler shift and often show no sig- 
nificant variation of the intensity and line width. 
The line profiles do not show any obvious asym- 
metry. These are probably spectroscopic signa- 
tures of kink/Alfven waves. In addition, we have 
presented a toy model, demonstrating that such 
transverse waves could also produce an observed 
intensity change and a tt/2 phase shift between 
the intensity and Doppler shift oscillations if the 
associated loops have monotonic density variation. 
Thus, the value of the phase shift alone is not a 
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reliable diagnostic of the wave mode in spectro- 
scopic observations made by instruments such as 
HINODE/EIS. 
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